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We present the main features of the spin-orbital superexchange which 
describes the magnetic and optical properties of Mott insulators with orbital 
degrees of freedom. In contrast to the SU(2) symmetry of spin superex- 
change, the orbital part of the superexchange obeys the lower cubic sym- 
metry of the lattice and is intrinsically frustrated. This intrinsic frustration 
and spin-orbital entanglement induce enhanced quantum fluctuations, and 
we point out a few situations where this leads to disordered states. Strong 
coupling between the spin and orbital degrees of freedom is discussed on 
the example of the RVOz perovskites, with R standing for rare-earth ion. 
La,- ■ ■,Lu. We explain the observed evolution of the orbital Too and Neel 
Tjvi transition temperature in the iiVOa series with decreasing ionic radius 
r_R. A few open problems and the current directions of research in the field 
of spin-orbital physics are pointed out. 
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1. Orbital versus spin superexchange 

In recent years the physical properties of Mott (or charge transfer) insu- 
lators are in the focus of interest in the condensed matter theory [1]. In order 
to develop theoretical understanding of complex phenomena in doped correlated 
insulators, including high temperature superconductivity in the cuprates and the 
colossal magnetoresistance (CMR) in the manganites, it is necessary to describe 
first the undoped materials, such as La2Cu04 and LaMnOs. In both cases the 
local Coulomb interactions (Hubbard U) is large and suppresses charge fluctua- 
tions, leading to low-energy effective Hamiltonians with superexchange interac- 
tions which stabilize antiferromagnetic (AF) spin order at low temperature [2, 3]. 
However, these two compounds are qualitatively quite different. On the one hand, 
the degeneracy of partly filled Cg orbitals is lifted in La2Cu04 by the tetragonal 
distortions of CuOg octahedra, resulting in two-dimensional (2D) AF superex- 
change of S' = 1/2 holes in — orbitals of Cu'^+ ions. On the other hand, in 
LaMnOa eg orbital degeneracy plays a fundamental role and, together with the 
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Jahn- Teller (JT) distortions, is required to understand the origin of the anisotropic 
^-type AF {A-AF) order [4]. As realized already three decades ago [2], in this 
latter case the orbital degrees of freedom, which are described by the components 

\a:)^l={x'-y^), |,)^_L(3.2_,2) (1) 

of pseudospin r = 1/2 and contribute explicitly to the structure of the superex- 
change. Thus they have to be included on equal footing with ionic spins in a 
spin-orbital superexchange model. In the last two decades several new concepts 
were developed, such as enhanced quantum fluctuations due to orbital degrees 
of freedom which participate in joint spin-orbital excitations [5], and spin-orbital 
entanglement which occurs in cases when spin and orbital operators cannot be 
decoupled from each other [6]. The actual physical problems in this emerging 
and rapidly developing field were reviewed in the Focus Issue Orbital Physics in 
New Journal of Physics [7] a few years ago. Here we want to focus on a few 
representative recent developments in the field of spin-orbital physics. 

Before we will discuss a few consequences of realistic spin-orbital models, let 
us consider first purely orbital superexchange interactions, as realized in ferromag- 
netic (FM) 2D systems, such as K2CUF4, or in FM ab planes of three-dimensional 
(3D) perovskitcs, either in KCuFa or in LaMnOa. The superexchange originates 
from charge d:-(fj ^ dfd^'^ excitations (with = t2ge^, df = t^g^l) the first 
case (for Cu^+ ions), and from d^dlj ^ dfd^ (with = t'lgCg, = t^g, and 
= ^ig'^g) excitations in the second one (for Mn'^^ ions) [8]. These charge exci- 
tations can be transformed away by applying a similar canonical transformation 
to the one leading to the t-J model [9], commonly used to describe the supercon- 
ducting cuprates. For the FM compound only high-spin excited states contribute, 
and in the strong-coupling regime [t <C U) one obtains the effective low-energy 
Hamiltonian describing orbital superexchange between pseudospins r = 1/2 due 
to Cg electrons (holes in KCuFs) in a cubic crystal [5, 10], 

Ho = Hj + H,. (2) 

the first term Hj in Eq. (2) describes superexchange interaction which depends 
on the superexchange constant 

[here t is the (dda) hopping element between two \z) orbitals along the c axis, 
and U is the intraorbital Coulomb interaction], and on Hund's exchange element 

,= ^. (4) 

This interaction parametrizes the multiplet structure of transition metal ions when 
anisotropy of Coulomb and exchange elements is neglected [8]. In the simplest 
case of d^ system, the singlet-triplet splitting between the two lowest excited 
states is 2r] [11]. 
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Following Rcf. [10], wc will consider below a 3D (perovskite) system with 
FM spin order and orbital interactions described by the Hamiltonian (2). The 
origin of intrinsic frustration in the orbital superexchange is best realized by 
considering its form — 

m ^ ^ 

where 7 = o, 6, c refers to the cubic axes, and r\ = 1/(1 — 3r]) follows from 
the energy of the high-spin charge excitation [8]. The pseudospin operators are 
defined as follows, 

rt''' = \{-a!±V3at), ^ = \<rt, (6) 

where cr^*'^'* are Pauli matrices and sign + (— ) is selected for a bond {ij) along a 

(6) axis. The pseudospin interactions oc rf'^Vj"' favor alternating orbitals (AOs) 
on each bond. They arc fimdamcntally different from the SU(2)-symmetric spin 
interactions, as they: (z) obey only lower cubic symmetry, [ii) are Ising-like, 
having only one component of the pseudospin interaction which favors pairs of 
orbitals oriented along the bond (z-likc) and in the; plane perpendicular to the 
bond (x-like), and {Hi) change their form when the cubic direction is changed. 
The interactions (5) look to be classical, but in fact due to the form of of the 
pseudospin Cg operators (6) they are not. However, the quantum corrections 
generated by them are rather small [10]. 

The second term in Eq. (2) stands for the orbital splitting (by Jsz), 
and follows from finite crystal field at either Cu^+ or Mn^+ ions, 

H, = -Je,Y.^'t\ (7) 

The crystal-field splitting favors electron (hole) occupancy in x [z] orbital for 
£2 > or (ez < 0). This term is induced by static distortions in the tetragonal 
field and is of particular importance in 2D systems, e.g. in the high- Tc cupratcs. 

The directional nature of pseudospin orbital interactions is responsible for 
their intrinsic frustration [5, 12]. In fact, the pair of orbitals which would minimize 
the energy for a bond (ij) is different for each cubic axis. Therefore, while for a 
single bond the minimal energy of —jJ can easily be obtained by selecting a pair 
of orthogonal orbitals on both sites, such as \z) and |a;) orbital for a bond along 
the c axis (note that this corresponds to the Ising superexchange interaction as 
only one orbital is active in intcrsitc charge excitations) , this is no longer possible 
for a 2D (or 3D) system. Thus, unlike in spin systems, the tendency towards 
disorder (quantum orbital liquid) is enhanced with increasing system dimension 
[13]. 

The essence of orbital frustration which characterizes the eg orbital superex- 
change (5) is captured by the 2D compass model, originally developed as a model 
for Mott insulators [2]. Intersite interactions in the compass model are descibed 



4 

by products t"t" of pseudospin components, where a = x,y, z, rather than by a 
scalar product Ti ■ Tj. 

(8) 

where a = x,y,z, rather than by a scalar product n ■ tj. In the 2D case the 
T^I^Tj interactions for bonds {ij) along the a axis compete with the t^t^ ones 
along the b axis [14]. Recently certain aspects of this model were investigated 
by analytic [14] and numerical [15, 16] methods. Despite its closeness to ordi- 
nary models used in quantum magnetism, there is no ordered phase with finite 
magnetization. Thus, this case is qualitatively different from the frustrated Ising 
interactions on the 2D square lattice, where it was proven rigorously that mag- 
netic order develops below certain transition temperature which depends on the 
ratio of frustrated to nonfrustrated plaquettes [17]. A competition of pseudospin 
interactions along different directions results instead in intersite correlations sim- 
ilar to the anisotropic XY model, and competition between two types of Ising-like 
order generates a quantum phase transition with high degeneracy of the ground 
state when all interactions have the same strength [15]. 

It is interesting to note that a similar quantum phase transition exists also 
in the one-dimensional (ID) compass model [18], when both above interactions 
alternate along the chain (A''' = A''/2 is the number of unit cells): 

N' 

n^D = U-^2i-ir^i + J^r^i^ii+i} ■ (9) 

This model was solved exactly in the whole range of {J^, Jz] parameters [18] by 
maping to the exactly solvable quantum Ising model [19] in different subspaces 
(see also this volume [20]). Equal coupling constants = Jz = J correspond to 
the the quantum critical point, where the orbital liquid emerges from two different 
disordered phases, with hidden order of pairs of pseudospins on every second bond. 
Presented here ID model model (9) provides thus a beautiful example of a first 
order quantum phase transition. 

2. Spin-orbital models and exotic liquid states 

The structure of spin-orbital superexchange involves interactions between 
SU(2)-symmetric spin scalar products Si ■ Sj on each bond (ij) connecting two 
nearest- neighbor transition metal ions, each one coupled to orbital operators 
{Ti,Tj} which obey only much lower symmetry (at most cubic for a cubic lat- 
tice) and its general form is [8], 

m 

The orbital operators J-j'^ and IC^'j'' involve the active orbitals on each bond 
(ij) along direction 7, either eg or t2g, which participate in dfd^ ^ d"~^^dj~'^ 
virtual excitations, and thus have the symmetry of the lattice {e.g. cubic in the 
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pcrovskitcs). The supcrcxchangc model (10) consists of several terms which orig- 
inate from different charge excitations. This feature made it possible to relate the 
averages of these different excitations to the spectral weights in the optical spec- 
troscopy [21], and this serves now as a theoretical tool to explain the observed 
anisotropy and temperature dependence of the spectral weights in the optical 
spectra [8]. 

The best known spin-orbital model is the Kugcl-Khomskii model for cP Cu^^ 
ions which involves the Cg orbital operators (6) and includes both FM and AF 
spin interactions, depending on whether the above charge excitation (for n = 9) 
involves the high-spin (triplet) or one of the low-spin (singlet) cP states. These 
excitations are parametrized by the coefficients ri = 1/(1 — 3r]), r2 = 1/(1 — 77) 
and r4 = 1/(1 + r/). Using compact notation of Ref. [11] it can be written as 
follows: 



-ri { S^-Sj + +r2 (s^-Sj - ^ 



1 __(7)_(7) 

4 ^ i 



m 

+(r.+r.)(l-ri^Vi^))(i-ri^))(i-.i^))}. (11) 

In ordeer to derive magnetic excitations for the systems with orbital de- 
generacy one usually derives magnetic exchange constants for a bond {ij) by 
averaging over the orbital operators in Eq. (10), 

Jij = {JS'')- (12) 

This procedure assumes implicitly that spin and orbital operators can be decou- 
pled from each other and ignores the possibility of entanglement and composite 
spin-orbital excitations introduced in Ref. [5]. It turns out that such excitations 
play a prominent role in destabilizing the classical AF long-range order in the 
spin-orbital model [22] . 

Qualitatively, the classical phase diagram of model (11) including the crystal- 
field term (7) consists in the {sz, if\ plane of four AF phases shown schematically 
in Fig. 1. First, if rj is small, then: (i) if £z > 0, the order is AF and the holes 
occupy \x) orbitals (2D AFx phase), while (ii) for < — \z) orbitals are oc- 
cupied in the anisotropic 3D AFz phase. Second, if Jh/U is large the high-spin 
excitations become important for the bonds with different orbitals on both sites, 
and weak FM interactions (oc r/J) occur within the ab planes, accompanied by 
strong (oc J) AF interactions along the c-axis, leading to ^-AF phase. Depending 
on the actual value of e^, either (Hi) A-AFl phase, or (iv) A-AF2 phase [with 
AF be planes and FM bonds along the a axis] is stable on the classical level. The 
energies of all the classical phases are degenerate at the (e^, 77) = (0, 0) point. The 
competition between the above types of spin and orbital order leads here to the 
quantum critical point at (e^,??) = (0,0). 

It has been argued [5, 22] that, similar to the situation in frustrated J1-J2- 
J3 spin models in quantum magnetism where the quantum corrections obtained 
from the spin-wave theory diverge and suppress the long-range order near quan- 
tum critical points [23], the spin order in all the G-AF and A-AF phases is sup- 
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Fig. 1. Schematic representation of orbital and magnetic long-range order in the ab 
plane of the (f model (11) in: G-AFx, G-AFz, yl-AFl, and A-AF2 phase. The shadded 
parts of different orbitals are oriented along the c-axis. The spins (arrows) in the next ab 
plane along the c axis arc AF to those below them in both G-AFz and A-AFl phase, and 
FM in the A-AF2 phase. In the G-AFx phase there is no superexchange coupling along 
the c axis to the next plane, but this degeneracy is easily removed when |a;) orbitals mix 
weakly with 1^;) orbitals. 

pressed near the (sz/J, 77) = (0,0) point. Although couterarguments were also 
given that the AFz (but not AFa;!) phase could be stabihzed by order-out-of- 
disorder mechanism [24] and this controversial issue awaits still an answer with 
a more sophisticated treatment, it is believed that only short-range order with 
correlations of valence-bond (VB) type may survive near the transitions between 
different phases with classical order (Fig. 1). In Ref. [5] a few representative 
variational VB wave functions (see Fig. 2) were constructed, with spin singlets 
on individual bonds {{ij)} accompanied by the variationally optimized orbitals 
pointing predominantly along the bond. Each VB phase is characterized by a 
set K which defines the bonds {{ij}} occupied by spin singlets, covering the 3D 
lattice. The two simplest (dimer) VB states are (see Fig. 2): (i) singlets along the 
a axis with orbitals close to Sx"^ — (VBa phase, degenerate with the analogous 
VB6 phase), (ii) singlets along the c-axis with \z) orbitals (VBc), are indeed more 
stable than the above classical AF phases, and we emphasize that this happens 
in three dimensions. 

Further improvement of the energy might be expected by including leading 
quantum fluctuations in the VB (VBa and VBc) states. In case of the VBa phase 
this would lead to resonance on each plaquette between the two components 
composed of spin singlets along the a and b axis, respectively, but this is hindered 
by the optimized z-likc orbitals, oriented in these states either along the a, or 
along the b axis (therefore the overlap between these two configurations is severly 
reduced). Unlike in the Heisenberg antiferromagnet, the bonds not occupied by 
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Fig. 2. Spin singlets (double lines) in VB states. The singlets along a axis (VBa) or 
along c axis (VBc) and ordered in either PVBA or PVBIc phase. The axes {a, 6, c} are 
shown by arrows. 

singlets contribute finite energy due to the orbital terms similar to those in Sec. I, 
and it turns out to be better to optimize the orbital energies due to those bonds 
depending on whether the bond connects two singlets: either (i) along a single 
(a, h, or c) axis, or {ii) oriented along two different axes with an angle of 7r/2, 
or {Hi) parallel to each other, with the bond making itself an angle of 7r/2 to 
both of them. As the second type of non-singlet bond is energetically the most 
favorable one, more energy is gained if plaquettes occupied by singlets along a and 
h axis alternate and form a 3D superlattice, a plaquette VB alternating (PVBA) 
state (Fig. 2), giving a lower energy that the quantum-corrected A-PiF phase. 
In contrast, the energy of the VBc phase is decreased more by the resonance 
between the vertical singlets {{ij) \\ c). The energy of the resulting resonating 
VBc (RVBc) state could be obtained using the Bethe ansatz result for the ID 
Heisenberg antiferromagnet, and adding the orbital energies due to the bonds 
{ij) II ah. 

Furthermore, a PVB interlayered phase along the c axis (PVBIc), composed 
of single planes of the PVBA phase interlayered with two planes of VBc vertical 
singlets (Fig. 2), is more stable in the crossover regime between the RVBc and 
PVBA phases for varepsilouz < 0. The orbital energies ex J in the bonds which 
connect the singlets on the bonds along the c axis with those lying in the ab planes 
are then optimized, while the energies cx are gained in the double layers of VBc 
phase. For varepsiloriz > 0, a rotated interlayered phase (PVBIa) is more stable 
before the PVBA takes over. Altogether, one finds [5] that an exotic spin-orbital 
liquid state, represented here by RVBc, PVBIc, PVBIa and PVBA phases, is 
favored in an extended regime of parameters, in analogy with a 2D 1/5-depleted 
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lattice [25], but the present instability is stronger and happens in a 3D modeL 

Finally, we remark that the present (f model (11) applies also to (f low-spin 
Ni^+ ions in LiNi02 via a particle-hole transformation. However, the situation is 
rather complex there and also other models were proposed. First, it was argued 
that a model based on symmetry arguments [26] , characterized by a large number 
of low-lying singlets associated to dimer coverings of the triangular lattice, could 
explain the properties of LiNi02. Second, it was shown that charge-transfer terms 
contribute to the superexchange and they change the balance between different 
terms in the Hamiltonian, making the orbital interactions stronger than the spin 
ones [27]. In any case, intcrplane JT interplanc coupling seems to be too weak 
in LiNi02 to stabilize the orbital long-range order. The observed difference in 
the physical properties between NaNi02 and LiNi02 remains one of the puzzling 
phenomena in the field [27] and awaits a future study which has to include all the 
above aspects. 

3. Entanglement in spin-orbital models 

Future developments in the theory to understand better the nature of or- 
dered and disordered states in spin-orbital systems have to take into account 
possible entanglement between both spin and orbital degrees of freedom on the 
bonds, as for instance in the states of Fig. 2. It is a common wisdom that 
the magnetism of correlated Mott insulators can be understood by means of the 
Goodenough-Kanamori rules. They originate from the mean-field (MF) picture 
and predict that the superexchange interaction between two magnetic ions with 
degenerate orbitals is strongly AF because of the Pauli principle if the overlap 
between the occupied orbitals is large, whereas it is weakly FM when the overlap 
is weak or virtually disappears. This means that spin order and orbital order 
are complementary — ferro orbital (FO) order supports AF order, while AO or- 
der supports FM spin order. Indeed, these celebrated rules are well followed in 
the manganites [28, 29], where the AO is robust in the FM ab planes, while the 
orbitals obey the FO order along the c axis, supporting the AF coupling and 
leading to the A-AF phase. In fact, one may use here the MF decoupling of the 
spin and orbital operators, which is sufficient to explain both the magnetic and 
optical properties of LaMnOa [8] . The orbital order (or liquid) state determines 
as well whether the intersite spin correlations are AF or FM in monolayer [30] 
and in bilayer [31] manganites, while in the perovskite Lai_a;Sra;Mn03 systems 
the magnons in the FM metallic phase are well explained by the orbital liquid 
state [32]. 

As a prominent example of the Goodenough-Kanamori complementarity we 
would like to mention also the AF phases realized in YVO3 [33], which are the 
subject of intense research in recent years. While xy orbitals of V'^^ ions are occu- 
pied due to GdFeOa distortions, the second electron of the cP ionic configuration 
may occupy either yz or zx orbital. Depending on whether these orbitals follow 
FO or AO order along the c axis, the magnetic correlations are there cither AF 
or FM, explaining the origin of the two observed types of AF order shown in Fig. 
3: (i) the C-AF phase, and (ii) the G-AF phase (with AF spin order along all 
three cubic directions). 
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Fig. 3. Goodenough-Kanamori rules on the example of two AF phases observed in the 
-RVO3 perovskites, with complementary magnetic and orbital order in the oc plane: (a) 
C-AF spin order accompanied by G-AO order; and (b) G-AF spin order accompanied by 
G-AO order. Arrows indicate spin order, while squares stand for two active tig orbitals, 
a and h. Both spins and orbitals alternate along the 6 axis (not shown). 

One may verify the Goodenough-Kanamori rules by evaluating intersite spin 
and orbital correlations, 

S,, = {S, ■ Sj)l{2Sf , T., = (f. . f,) , (13) 

and comparing them with each other. In addition, spin-orbital entanglement can 
be measured by the composite correlation function defined for a bond (ij) [6], 

= {{{S. ■ • T,)) ~ (5,; • S.^{% ■ f,)]/i2Sf . (14) 

When Cij = 0, the spin and orbital operators are disentangled and their MF 
decoupling can be applied, while if Cij < — spin and orbital operators are 
entangled. Consider again the structure of spin-orbital superexchange models 
(10) derived for t2g systems from charge excitations between transition metal ions 
in configurations: d} (Ti^+ ions in cubic titanates) with S = 1/2 [34], and (f {Y^+ 
ions in cubic vanadates) with S = 1 [35], for more details about the structure of 
Hj in both and model see for instance Ref. [8]). 

To capture the essence of spin-orbital entanglement it is sufficient to solve 
t2g models for d^ and d^ configurations on four-site chains oriented along the c 
axis using periodic boundary conditions. One finds entangled spin-orbital states 
in both models, and nontrivial spin-orbital dynamics strongly influences both 
intersite spin and orbital intersite correlations (13). In the d^ (titanate) case 
one recovers the SU(4) model with Sij = = dj = —0.25 at 77 = [Fig. 
4(a)]. At finite rj the SU(4) degeneracy of all intersite correlations is removed and 
Tij < dj < Sij < 0, as long as the spin singlet {S = 0) ground state persists, i.e., 
for < 0.21. In this regime the Goodenough-Kanamori rule with complementary 
correlations (implying different signs of Sij and Tij), is violated. Instead, by 
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Fig. 4. Left — intersite spin Sij = {Si-Sj) (filled circles), orbital Tij = (Ti-Tj) (empty 
circles), and composite spin-orbital dj (f4) (crosses) correlations for increasing Hund's 
exchange r]S along c axis; — right — the corresponding spin exchange constants Jij (14), 
as obtained for A'^ = 4 site spin-orbital chain (10) with PBC for: (a),(c) (titanate, 
S = 1/2) model, and (b),(d) d? (vanadate, S = 1) model. In the shaded areas of (c) 
and (d) the spin correlations Sij < do not follow the sign of the exchange constant 
Jij < 0, and the classical Goodenough-Kanamori rule is violated. 



analyzing the values of the intersite correlations one finds that the ground state 
wave function for each bond (ij) is close to a total spin-orbital singlet, involving a 
linear combination of (spin singlet/orbital triplet) and (spin triplet/orbital singlet) 
states. The vanadate (P model behaves in a similar way, with all three Sij, Tij and 
Cij correlations being negative in the spin-singlet (S = 0) regime of fluctuating 
yz and zx orbitals obtained for rj < 0.07 [Fig. 4(b)]. Therefore, the composite 
spin-orbital correlations are here finite {Cij < 0), spin and orbital variables are 
entangled, and the MF factorization of the ground state fails. 

To provide further evidence that the Goodenough-Kanamori rules do not 
apply to t2g systems in the regime of small rj, we compare spin exchange constants 
Jij (12) with the actual values of intersite spin correlations Sij (13). One finds that 
exchange interaction is formally FM {Jij < 0) in the orbital-disordered phase at 
low values of rj [Figs. 4(c) and 4(d)] is in fact accompanied by AF spin correlations 
{Sij < 0), so Jij Sij > and the ground state energy calculated in the MF theory is 
enhanced [6]. In contrast, similar analysis (not shown) performed for spin-orbital 
model (11), derived for ions with eg orbital degrees of freedom, gave Jij Sij < 
0, so here spin correlations follow the sign of the exchange constant [6]. This 
remarkable difference between t2g and Sg systems originates from composite spin- 
orbital fluctuations, which are responsible for the 'dynamical' nature of exchange 
constants in the former case. They exhibit large fluctuations around the average 
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value, measured by 5J = {{{jlj^'Y) — Jfj}^^"^- As an illustrative example, we 
give here the values found in the and cP model at 77 = 0. While the average 
spin exchange constant is small in both cases {Jij ~ for d^, Jij ~ —0.04 for rf^), 
jj^j'^ fluctuates widely over both positive and negative values. In the d^ case the 
fluctuations between {S = 0/T — 1) and {S = 1/T = 0) bond wave functions are 
so large that SJ = l\ They survive even quite far from the high-symmetry SU(4) 
point (at T] > 0.1), and stabilize spin-orbital singlet phase in a broad regime of 
T]. Also in the d^ model the orbital bond correlations change dynamically from 
singlet to triplet, resulting in dJ > | |, with SJ = i{l - (2T,j + \fYl'^ ~ 0.25, 
while these fluctuations are small for model (11) with eg orbitals. 

Wc emphasize that composite spin-orbital fluctuations which occur in spin- 
orbital entangled states for realistic parameters determine the magnetic and op- 
tical properties of titanates and vanadates. For instance, such composite spin- 
orbital fluctuations are responsible for the temperature dependence of the optical 
spectral weights in LaVOs [21] and trigger spin-orbital dimerization in the C-AF 
phase of YVO3 in the intermediate temperature regime (see below) [36]. Re- 
markably, the observed dimerization in the magnetic excitations may be sees as a 
signature of entanglement in excited states which become relevant at finite tem- 
perature. It is activated by thermal fluctuations in the spin chain [37], which 
couple to dimerized correlations in the orbital sector. 

4. Phase transitions in the iJVOs perovskites 

Recent progress in experimental studies of transition metal oxides provided 

exceptionally detailed information on the phase diagrams of the iJMnOs mangan- 
ites [38] and the RVO3 vanadates [39] (where i?=La,Pr,- • •,Lu). In the manganites 
the orbital order appears first when the temperature is lowered (at Too ^ 800 
K), and its onset is accompanied by a lattice distortion. The magnetic order 
follows at much lower temperature Tn ^ 140 K. This separation of the spin and 
orbital energy scales follows also from the JT distortions which contribute to the 
orbital interactions and enhance the value of Tqq [4]. Recent experiments showed 
that the orbital transition temperature Too is enhanced when the ionic radius rn 
of the ions decreases along the i?Mn03 perovskites, while is drastically 
reduced, resulting in the change of magnetic order from the j4-AF to the ii^-AF 
phase [38]. This behavior could not be understood until now and is one of the 
challenges for the theory. 

The expcrimcintal situation in the cubic vanadates is different and even more 
complex. On the one hand, the magnetic order in YVO3 is G-AF [Fig. 5(b)] at 
low temperature, and changes at the first order magnetic transition at = 77 
K to the C-AF structure, which remains stable up to T^vi ~ 116 K. The magnetic 
transition at T^2 is particularly surprising as the staggered moments change their 
direction from approximately parallel to the c axis in the G-AF phase to lying 
almost within the ab planes in the G-AF phase, with some small alternating G- 
AF component [40]. Although it was argued that the entropy due to magnetic 
and orbital excitations is higher in the G-AF phase [35, 41], the magnetization 
reversal at the lower magnetic transition remains misterious, and this phenomenon 
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is still puzzling and far from being completely understood. On the other hand, 
only the C-AF order develops in LaVOa below a somewhat higher T/vi — 143 
K, and is almost immediately followed by a weak structural transition stabilizing 
the weak G-AO order at Too - 141 K [39, 42] [Fig. 5(a)]. Remarkably, the 
magnetic order parameter in the C-AF phase of LaVOa is strongly reduced to 
~ 1.3/iB, much below the reduction expected from quantum fluctuations in the 
C-AF phase (being only 6% for S' = 1 spins) [43] — also this reduction of the 
measured magnetization could not be explained so far. 

Experimental studies have shown that the C-AF order is common to the en- 
tire family of the i?V03 vanadates, and in general T/vi < Too- except for LaVOa 
with Tni ~ Too [39, 42]. When the ionic radius vr decreases, the Neel tem- 
perature Tjvi also decreases, while orbital transition temperature Too increases, 
passes through a maximum close to YVO3, and decreases again towards LuVOa. 
Knowing that orbital quantum fluctuations and spin-orbital entanglement play 
so important role in the perovskite vanadates, it is of interest to ask whether 
the spin-orbital model of the form (10) introduced for the perovskite vanadates 
[35, 36] is able to describe this variation of Too and Tati with decreasing radius 
rn of R ions in KVO3 [39] . It is clear that the nonmonotonic dependence of Too 
on rji cannot be reproduced just by the superexchange, as it requires two mech- 
anisms which oppose each other. In fact, there is even no reason to assume that 
the superexchange constant J should depend on r^, as the distances between 
ions are quite close to each other in different compounds [44, 45, 46], so neither t 
nor U is expected to change significantly. 

As in the iJMnOs manganites [38], one expects that the JT distortions 
should increase when the ionic radius r/? decreases, and one could argue that this 
would induce the increase of Too- In order to unravel the physical mechanism 
responsible for the decrease of Too from YVO3 to LuVOs one has to analyze in 
more detail the evolution of GdFeOs distortions with decreasing rji [47]. Such 
distortions are common for the perovskites [48], and may be described by two 
subsequent rotations of VOe octahedra: (i) by an angle 1? around the b axis, and 
(ii) by an angle ip around the c axis. Increasing angle d causes a decrease of V O- 
V bond angles along the c direction, being tt — 2'd, and leads to an orthorhombic 
lattice distortion u = (b — a) /a, where a and b are the lattice parameters of the 
Pbnm structure of i?V03 . The structural data for the perovskite i?V03 vanadates 
[44, 45, 46] give the following empirical relation between the ionic radius tr and 
the angle 1?: 

rR = ro-asm'^'&, (15) 

where rg = 1.5 A and a = 0.95 A. This allows one to use the angle '& to parametrize 
the dependence of the microscopic parameters and the transition temperatures 
Too and Tni on vr. 

The spin-orbital model introduced in Ref. [47] to describe the phase dia- 
gram of i?V03 includes: (i) the superexchange cx J between V^^ ions in the (P 
configuration with 5 = 1 spins [35], (ii) the crystal-field splitting cx between 
yz and zx orbitals, (iii) intersite orbital interactions oc Vab,Vc (which originate 
from the coupling to the lattice), and {iv) orbital-lattice term which induces or- 
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bital polarization when the distortion u inreases. The Hamiltonian consists thus 

of several terms [47], 

(«J>lt7 « {ij)\\c 

+ Vab{^)Y,r:r^^-guY,rr + lNK{u-uom^. (16) 

{ij)\\ab i 

with 7 = a, 6, c labels the cubic axes. The orbital operators take the form: 
Jij^ = ^{(l + 2r?ri)(^fi-f, + in,n^.^ 

- Vrs (^n X fj + ^n.n^^ - ^imim + nj)| , (17) 
K-]^ = (^n ■ fj + ^riiTi^ + r]r3 (ji x fj + in^n^- 



1 



(7) 



-{l + rjri){ni + nj)j , (18) 

and have a rich structure which originates from the projections of the excited 
states on the respective eigenstates of ion. They arise due to the ^ dfdj 
charge excitations, leading either to high-spin or to low-spin configurations, so 
Hund's exchange splittings in the multiplet structure of a V^^ ions enter via 
the coefficients ri = 1/(1 — 3r]) and rs = 1/(1 + 2??) (the lower energy singlet 
excitations occur at energy ?7, so the corresponding coefficient is r2 = 1). Using 
the atomic value of Hund's exchange, we have estimated that j] ~ 0.13 [8, 41]. The 
orbital operators fi = {t^ ,t~ ,tI} for pseudospin r = 1/2 in Eqs. (17) and (18) 
are defined in the subspace spanned by two orbital flavors which are active along 
a given cubic direction 7 = a,b,c. Unlike in (P model (11) for Cg orbitals, the 
leading orbital interactions are here proportional to the scalar products {fi ■ fj)^^"^ 
of orbital operators on the bonds, but the structure of local Coulomb interactions 
is responsible for additional terms, 

in X fj)^^^ = I (r+r/ + r-r-) + r^r^ , (19) 

so the orbital quantum numbers are not conserved. Finally, the operator n^^-* 
stands for the number of active electrons at site i along the bond (ij) , for instance 
for a bond along the c axis this number is n^f^ = ^iyz ^ zx ■ 

Further insight into the electronic configuration on V'^^ ions in i?V03 may 
be obtained either by investigating the electronic structure [49], or by explicit 
calculations using the point charge model [47]. Both lead to the electronic con- 
figurations {xyY{yz / zxY , i.e., the xy orbitals are occupied at all sites, while the 
remaining {yz^ zx} orbitals represent active orbital degrees of freedom which con- 
tribute to the t2g orbital dynamics, expressed by the scalar product f ■ fj (only 
for 7 = c). This demonstrates an important difference between the eg (with one 
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electron or one hole in Cg orbitals at each site) and a t2g system, such as the i?V03 
perovskites. Two active t2g orbitals along this bond open a new possibility — if 
both orbitals are singly occupied, an orbital singlet gives here FM superexchange, 
even in the absence of Hund's exchange (at = 0) [35]. 

Furthermore, the actual electronic configurations realized in the KVO3 sys- 
tems justify the form of the crystal-field term used in Eq. (16), with r? = 
^{^i^yz — ni^zx)- The crystal field alternates in the ab planes, but is uniform 
along the c axis. It is thus characterized by the vector Q = (tt, tt, 0) in reciprocal 
space, and competes with the (weak) G-AO order supporting the observed C-AF 
phase below Tni- The orbital interactions induced by the distortions of the VOe 
octahedra and by GdFeOs distortions of the lattice, Vab > and 14 > 0, also 
favor the C-AO order (like > 0). Note that 14 > counteracts the orbital 
interactions included in the superexchange. In case of K one may assume that its 
dependence on the angle is weak, and a constant T4(i?) = 0.26 J was chosen in 
Ref. [47] in order to satisfy the experimental condition that the C-AF and G-AO 
order appears almost simultaneously in LaVOs, sec Fig. 5(a). The experimental 
value = 1^3 K for LaVOs [39] was fairly well reproduced in the present model 
taking J = 200 K. The functional dependence of the remaining two parameters 
{Ez,Vab} on the tilting angle d was derived from the point charge model [47], 
using the structural data for RYO3 [44, 45, 46]. In this way the following two 
relations were deduced: 



The last two terms in Eq. (16) describe the linear coupling a g' > between 
active {yz, zx} orbitals and the orthorhombic lattice distortion u. The elastic 
energy which counteracts lattice distortion u is given the force constant K, and 
N is the number of V^~^ ions. The coupling oc gu acts as a transverse field in the 
pseudospin space. While the eigenstates -^(jy-^) ± \zx)) favored by rf cannot 
be realized due to the competition with all the other terms, increasing lattice 
distortion u (increasing angle -d) modifies the orbital order and intersite orbital 
correlations. At the minimum one finds an effective coupling, 



with the global distortion w(i?; T) caused by («) a pure lattice term Mo(t^)) and {ii) 
an orbital contribution oc (r^) which can be determined self-consistently within 
the present approach, see Ref. [47]. Both uq and (r^) are expected to increase 
with increasing tilting ■&. This dependence has to be faster than the ones of Eqs. 
(20) and (21), as otherwise no maximum in the dependence of Too on m would 
be obtained. Using this argument a semiempirical relation. 



E^{i^) = Jt;^ sin^i?cosi?, 
Vabi-^) = Jvab sin^'d cos -d. 



(20) 
(21) 



(22) 



Qeffi'^) = Jvg sin' 



1? cos 1? , 



(23) 



was postulated in Ref. [47]. Therefore, magnetic and orbital correlations de- 
scribed by spin-orbital model (10) depend on three parameters: {vz,Vab,Vg}. 
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Fig. 5. Theoretical description of orbital and magnetic phase transition in the RVO3 
perovskites obtained using model (16), for varying ionic radius vr of R ions: (a) tran- 
sition temperatures — Too for the onset of orbital and Tmi for the onset of the C-AF 
order, as obtained in Ref. [47] (solid lines), and compared with the experimental data 
of Ref. [39] (filled diamonds and open circles); (b) evolution of orbital order parameter 
(rf )g (filled circles), transverse orbital polarization (rf ) (squares), and orbital intersite 
correlations |(Ti-7^ )| (diamonds) along c axis at T = 0. Parameters: J = 200 K, Vz = 17, 
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22, 



740. 



While the above fast dependence on the tilting angle of VOe octahedra 
in the i?V03 family was introduced in order to reproduce the experimentally 
observed dependence of Too on r/j, see Fig. 5(a), it may be justified a posteriori. 
It turns out that the dependence of ^eff on the ionic radius r/j follows the actual 
lattice distortion u in i?V03 measured at T = (or just above T^i) [47]. Also the 
orbital polarization (r^) is approximately cx sin^ 1} cos 1}, and follows the same fast 
dependence of 5eff(i^) for the RVO3 perovskites. These results indicate that the 
bare coupling parameters {g, K} are nearly constant and independent of r^, which 
may be treated as a prediction of the theory to be verified by future experiments. 

Similar to the i?Mn03 manganites [4], the correct MF treatment of the 
orbital and magnetic phase transitions in the i?V03 vanadates requires the cou- 
pling between the on-site orbital, {t^)g = — ■^/)l> and spin, {Sf)c, order 
parameters in the C-AF phase, as well as a composite {S^rf) order parameter. 
Unlike in the perovskite manganites [38], however, the on-site MF theory does not 
suffice for the i?V03 compounds as the orbital singlet correlations {ri ■ fj) on the 
bonds {ij) along the c axis play an important role in stabilizing the C-AF phase 
[35, 41]. Therefore, the minimal physically acceptable approach to the present 
problem is a self-consistent calculation for a bond {ij) along the c axis, embedded 
in the MF terms due to its neighbors along all the cubic axes [47] . This procedure, 
with properly selected model parameters, served to reproduce the experimental 
phase diagram of Fig. 5(a). One finds that the orbital transition occurs first in 
the i?V03 vanadates with < rLaj as observed. At the same time, the Neel 
temperature T/vi decreases with decreasing r^. 
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The key feature of the present spin-orbital system which drives the observed 
dependence of Too and Tjvi on the ionoc radius tr is the evolution of intersite 
orbital correlations {ri-fj) along the c axis, and the increasing orbital polarization 
(r^) with decreasing tr [Fig. 5(b)]. Indeed, the singlet correlations are drastically 
suppressed from LaVOa towards LuVOa. While {t^) ~ 0.03 is rather weak in 
LaVOa, it steadily increases along the iJVOs perovskites when m decreases, and 
finally it becomes as important as the orbital order parameter, (r^) ^ {t^)g- 
Note that in all the cases the latter order parameter is substantially reduced from 
the classical value ^ by singlet orbital fluctuations, being {t^)g — 0.32 and 0.36 
for LaVOs and LuVOs, respectively. 

The results presented in Ref. [47] demonstrate the remarkable dependence 
of both spin and orbital order on the orbital correlations. First, the increase of 
orbital intersite interactions due to the JT term (21), induces steady increase of 
the orbital temperature Too with decreasing tr. When the orbital polarization 
(r^) becomes large, however, this increase is suppressed and one reproduces the 
nonmonotonic dependence of Too on tr, with the observed drop of Too when tr 
decreases beyond tr ~ 1.18 A of YVO3. Second, the changes in intersite orbital 
correlations shown in Fig. 5(b) modify the magnetic exchange constants {Jab, Jc} 
along the bonds parallel to one of cubic directions in the ah planes or along the 
c axis (12), and thus the value of Tjvi is reduced with decreasing tr, although 
the superexchange energy J does not change. This also means that the width of 
the magnon band given at T = by M^c-AF = '^{Jab + \Jc\) [41] is reduced by 
a factor close to 1.8 from LaVOa to YVO3, in agreement with surprisingly low 
magnon energies observed in the C-AF phase of YVO3 [33] . 

Summarizing, spin-orbital model (16) of Ref. [47] provides an almost quan- 
titative understanding of the systematic experimental trends for both orbital and 
magnetic transitions in the iJVOs perovskites, and is able to reproduce the non- 
monotonic variation of the orbital temperature Too for decreasing tr. Hovever, 
the theoretical description of the magnetic transition to the G-AF phase at Tjv2, 
which occurs for small tr [39], remains to be addressed by future theory. More 
open issues and future directions of reasearch in the field of the orbital physics 
are shortly indicated in the next section. 

5. Summary and open problems 

A few representative problems discussed above demonstrate that the spin- 
orbital physics is a very rich field, with intrinsically frustatred interactions and 
rather exotic ordered or disordered phases, and with their behavior frequently 
dominated by quantum fluctuations. While valuable information about the elec- 
tronic structure is obtained from density functional theory [50], the many-body 
aspects have to be studied using models of correlated electrons. Perhaps the most 
important feature of orbital or spin-orbital superexchange models is strong cou- 
pling between the orbitals and the lattice, which may help to stabilize orbital 
order, as in model (16) introduced for the i?V03 perovskites [47]. In this way, 
the lattice distortions may also indirectly influence the onset of magnetic order 
in systems with active orbital degrees of freedom. If they are absent and the 
lattice is frustrated in addition, a very interesting situation arises, with strong 
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tendency towards truly exotic quantum states [51]. Examples of this behavior 
were considered recently for the triangular lattice, both for eg orbitals in LiNi02 
[27] and t2g orbitals in NaTi02 [52] . None of these models could really be solved, 
but generic tendency towards dimcr correlations with spin singlets on the bonds 
for particular orbital states could be shown. 

Rapid progress is the field of orbital physics results mainly from developing 
new experimental techniques and synthesizing novel materials. While the exper- 
iment is ahead of the theory in most cases, there are some exceptions. One of 
them was a theoretical prediction of the energy and dispersion of orbital excita- 
tions [10, 53]. Only recently orbital excitations (orbitons) could be observed in 
Raman scattering in the Mott insulators LaTiOa and YTiOa [54] . An interesting 
question which arises in this context is carrier propagation in a Mott insulator 
with AO order. This problem is very complex and has been addressed also in the 
last decade of the previous century [55]. Not only spin excitations may provide 
more final states, as for triplet excitations in the t-J model [56], but a carrier 
may dress by orbital excitations [57], and scatter on magnons [58]. Indeed, the 
coupling to orbitons increases the effective mass of a moving hole in Cg systems 
[57]. The orbital part of the superexchange is classical t2g systems, but also there 
weak quasiparticle dispersion arises from three-site processes [59] (see also this 
volume [60]). 

Similar to doped manganites [61], also in doped i?i_x(Sr,Ca)a:V03 systems 
orbital order gradually disappears [62]. The composite C-AF/G-AO order sur- 
vives, however, is a broad range of doping, in contrast to Lai-a^Sr^^MnOa, where 
FM order sets in already at a; 0.10. It is quite remarkable that the comple- 
mentary G-AF/C-AO type of order is fragile and disappears in Yi-ajCanVOa 
already at a; = 0.02 [62]. We remark that the hole doped in C-AF/G-AO phase 
are localized in polaronlike states [63], so the purely electronic model such as 
of Ref. [59] is too crude to capture both the evolution of the spin-orbital order 
in doped vanadates and the gradual decrease of the energy scale for spin-orbital 
fiuctuations. 

As a final remark, we would like to mention very promissing recent experi- 
mental studies of Ni-based superlattices [64]. Recent theory for LaNiOa/LaMOa 
superlattices (with M=A1, Gd, Ti) [65] predicts that the correlated eg electrons 
in the Ni02 planes develop a planar — orbital order driven by the reduced 
dimensionality and further supported by epitaxial strain from the substrate. This 
resembles the — orbital polarization in doped layered manganites [66] . As 
in all other cases discussed above, the superexchange interactions which involve 
the orbital degrees of freedom play here a crucial role to understand the observed 
magnetic order and low-lying excited states. 
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